Adatoms in graphene as a source of current polarization: Role of the local magnetic 
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We theoretically investigate spin-resolved currents flowing in large-area graphene, with and with- 
out defects, doped with single atoms of noble metals {Cu, Ag and Au) and 3d-transition metals 
{Mn,Fe,Co and Ni). We show that the presence of a local magnetic moment is a necessary but not 
sufficient condition to have a non zero current polarization. An essential requirement is the pres- 
ence of spin-split localized levels near the Fermi energy that strongly hybridize with the graphene tt 
bands. We also show that a gate potential can be used to tune the energy of these localized levels, 
leading to an external way to control the degree of spin-polarized current without the application 
of a magnetic field. 

PACS numbers: 72.80.Vp,73.23.Ad,71.15.Mb 
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I. INTRODUCTION 

Grapheneii^ possesses a very long mean free path of 
the order of micrometers and a ballistic transport regime 
at low temperatures^, which combined with a low spin- 
orbit coupling make it ideal for applications in spintronic 
devices. However, due to the intrinsic spin-unpolarized 
electronic structure of a pristine sheet, it is necessary 
to modify the electronic structure of graphene in such a 
way that the transport properties of the two spin chan- 
nels become different. Few ways that have been proposed 
to do this arc (i) quantum confinement (e. g. graphene 
nanoribbons^"— , nanoflakes^, hidrogenation^, and inter- 
action with substratesiSiii) ; ii) straini^; and iii) adsorp- 
tion of molecules or atom o^^d"^ , where the presence of a 
local magnetic moment leads to a spin-non-degenerated 
electronic structure. However, edge disorder is a limit- 
ing factor to use quantum confinement as a way to dif- 
ferentiate the spin channels, and for the application of 
strain, there is the drawback of a mechanically control- 
lable device. Therefore, doping large-area graphene with 
magnetic dopants seems more suitable for applications in 
spintronics. 

In the last few years, motivated by spintronic applica- 
tions, many studies focused on understanding the dop- 
ing properties of metallic adatoms on graphene. Experi- 
mentally, it has been shown that both noble metals and 
transition metals tend to diffuse and clusterizc in pristine 
graphene, unless they arc trapped in lattice defects such 
as monovacancies (MVs) and divacancies (DVs)J^"— 
Moreover, all theoretical studies of metallic doping in 
large-area graphene have focused solely in investigating 
the emergence of local magnetic momentsi^Ti^^. However, 
conclusions about the possibility of using metallic doping 
to generate spin-polarized transport requires the investi- 
gation of spin-resolved currents, instead of merely local 
magnetic moments. 

In this paper, we theoretically investigate spin-resolved 
currents fiowing in large-area graphene, doped with single 




FIG. 1: (Color Online) (a) Representation of devices consid- 
ered in this work, (b) Possible adsorption sites for atoms on 
pristine graphene. (c) Geometry for Co@MV, similar to all 
other atoms, (d) Co@585, similar to all other atoms, (e) 
Co@555777 (similar to Mn:Co:Ni@555777). (f) Au@555777. 
(g) Cu@555777 (similar to ^5@555777). 



atoms of noble metals (Cu, Ag and Au) and Sd-transition 
metals {Mn,Fe,Co and Ni). We study the adsorption 
both in pristine graphene as well as in monovacancies 
and divacancies, and we show that the current polariza- 
tion is not monotonically correlated to the local magnetic 
moment. Instead, it depends on the presence of localized 
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levels near the Fermi energy which strongly couple (and 
thus hybridize) to the graphene tt bands. We also show 
that a gate potential can be used to tune the energy 
of these localized levels, leading to an external way to 
control the degree of spin-polarized current without the 
application of a magnetic field. 

In Fig. [T] (a) we show the setup of the devices we 
considered in this work, and in panels (b)-(g) we show 
the representative geometries for the complexes of met- 
als adsorbed in graphene, both pristine and with defects 
(metal@defect). The paper is organized as follows: in the 
next section, we describe the methodology we used. In 
the subsequent section, we present our results, discussing 
the total energy calculations, the current polarization, 
the dependence on the doping concentration, and the ef- 
fect of a gate potential. In the last section we present 
our conclusions. 



II. METHODOLOGY 
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FIG. 2: (Color Online) Schematic representation of the uni- 
tary cell in the real (a) and reciprocal (b) spaces, where we 
define the symmetry points used to plot the energy bands. 



In order to obtain the geometries and electronic struc- 
ture, we employ total energy density functional theory^ 
(DFT) calculations. All geometries were fully relaxed 
using a force criteria of 0.02 eV/A with the Siesta 
code^. All calculations were carried out applying the 
generalized gradient approximation as parameterized by 
Pcrdcw, Burke, and Ernzerho^ (GGA-PBE) for the ex- 
change correlation functional (Exc)- We do not consider 



the spin-orbit coupling, because such interaction may lo- 
cally open small band gaps^^i^^, of the order of few milli 
electron volts, which will be a small effect when compared 
to the single-impurity scattering effects investigated in 
this paper. A localized double-C polarized basis, an elec- 
tronic temperature of bmeV, and 300 Ry for the grid 
cut-off were used. Most of our results were obtained with 
a rectangular 7x6 graphene supcrcell, as schematically 
shown in Fig. [2ja). The integration in the Brillouin zone 
was done with the Monkhorst-Pack^ scheme correspond- 
ing to a 48 X 48 sampling in the unitary cell. A vacuum of 
20A was used between the periodic images of graphene 
sheets to avoid undesirable interactions. The directions 
between the symmetry points defined in Fig. [5] (b) were 
used to plot the energy bands. 

Once the relaxed geometries were obtained, we cal- 
culated the transport properties with the Transampa 
code^. This code makes use of the non equilib- 
rium Green's function techniques coupled to DFT in 
a fully self-consistent way with the Biittiker-Landauer 
formalisni^ within the non interacting approach of Meir 
Wingreen^. Within this formalism, the translational 
symmetry along the transport direction is broken. How- 
ever, in the transversal direction, due to the period- 
icity of the infinite 2D-system, the Bloch theorem is 
still usable. In order to take into account a truly 2D 
graphene sheet, we implement the transport equations 
with transversal periodicity, considering fc-points perpen- 
dicular to the transport direction (fc_L). In the transport 
calculations, the relaxed geometries were sandwiched be- 
tween two buffer layers to guarantee that the scattering 
region has a pristine-like coupling with the leads. We 
used the same simulation parameters of the Siesta cal- 
culations, with exception of the /c-points sampling. Here, 
for the electronic density, we used the equivalent of 48 
/cj_-points in the unitary cell, whereas for the transmit- 
tance. 5000 /c^-points are necessary for a smooth behav- 
ior. This number of A:-points is sufficient to capture the 
true 2D-character of graphene. In each step of the self- 
consistent field cycle, we calculated the electronic density 
with an integral in the complex plane discretized in 60 
energy points and 5 poles. When a bias potential was ap- 
plied, the non equilibrium contribution for the electronic 
density was calculated with an integral discretized in 200 
energy points. 

We also implement the effect of a gate potential with 
strength "Vgate" by adding to the Hamiltonian a function 
V{z) that depends only on the z coordinate, given by: 







V{z) 
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Vgate 
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, z < zi or z > Z2 
I ,0<z-zi<r 

, Zl+T<Z<Z2-T 

> , — r < z — Z2 < 0. 



(1) 
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Since z is the transport direction, the above function af- 
fects only a region between z\ and Z2 (gated region) and 
is continuous up to its first derivative. This region en- 
closes the metal and the defect in our simulations, and 
we verified that the buffer-layer regions were not affected 
by this gate. The r parameter controls the smoothness 
of the connections between the gated and non-gated re- 
gions, and we used r = 1.25 A for a smooth behavior. 
It is also important to mention that the gate effect is in- 
cluded in the whole self-consistcnt-ficld procedure, and 
the electronic density rearrangement creates a screening 



effect, in the sense that (for example) 



gate 
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not necessarily shift the energy levels by leV . This im- 
plementation is in the same spirit of Refs. IH andH. 



III. RESULTS AND DISCUSSION 

A. Geometry, magnetic moment, and binding 
energies 

The relaxed geometries, local magnetic moments (m) 
and binding energies (i^h) are shown in table HI Our re- 
sults compare well with previous calculationsi^i^Si"— i^. 
This is also true for the adsorption of noble metals on 
pristine graphene, where van der Waals interactions arc 
important^. However, in these cases, the energies ob- 
tained for the different configurations with our method- 
ology arc all very similar. Thus, to select which adsorp- 
tion site to use in the charge transport calculations [B 
position in Fig. [Kb)], we used the results of Ref. [s^- For 
atoms trapped in the MV, all relaxations lead to geome- 
tries similar to Co@MV [depicted in Fig. [T] (c)]. How- 
ever, the transition metals have a C31, local symmetry, 
whereas for the noble metals this symmetry is destroyed 
by a structural distortion, breaking the spin degenera- 
cies of localized states near the Fermi level, leading to 
m = 1.0 p,B, as discussed in Ref. [2l|. For adsorption on 
DVs, we considered two particular reconstructions stud- 
ied in the literature^, the 585 and 555777. For atoms 
trapped in 585 DVs, we find a C2V local symmetry in all 
cases, with the atoms displaced from the graphene plane, 
similar to Co@585, as shown in Fig. [T] (d). All the id 
transition metals adsorbed on the 555777 DV present a 
non-zero local magnetic moment, with the energetically 
most favorable position being above the central atom of 
this defect, as shown in Fig. [T] (e). It is worth noting 
that (i) Ni has a non zero local magnetic moment solely 
when trapped in this site and (ii) adsorption of adatoms 
has been experimentally observed in thcniii. The noble 
metals, on the other hand, have a completely different 
behavior, presenting a zero local magnetic moment when 
adsorbed on the 555777 DV. The lowest energy config- 
uration for Ait@555777 is shown in Fig. [1] (f), with the 
gold atom at the top of a carbon atom belonging to a 
pentagon. For Cu and Ag^ the lowest energy position is 
above the center of a heptagon, as shown in Fig. [1] (g). 
However, Cu presents a Cg local symmetry, whereas Ag 
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FIG. 3: (Color online) Current polarization with Vbias = 
50meV (a) for all studied cases, (b) as a function of the local 
magnetic moment, (c) as a function of the magnetic moment 
induced at the graphene carbon atoms. The solid lines guide 
the eye. 



presents five different chemical bond lengths. 



B. Current Polarization 

In Fig. [3] (a) we show the current polariza- 
tion, CP=(/^ — /|)/(/f -I- /j,), for ah studied cases, 
where the current for each spin channel is I a = 
[e/h) J T„{E)[f{E - fiL) - f{E - fiB)]dE, with T,iE) 
the transmittance for the spin cr, f{E) the Fermi-Dirac 
function, and fJ.L{R) the chemical potential for the left 
(right) electrodes. A bias potential of 50 meV was used. 
Clearly, Co and Fe adsorbed on pristine graphene gen- 
erate the higher CP in the systems we consider. How- 
ever, the metal diffusion experimentally observed creates 
a drawback to build devices with adatoms over pristine 
graphene since in the most cases the atoms tend to clus- 
terize. For the atoms trapped in defects, the higher val- 
ues of CP (in modulus) occurs for the 7Vz@555777, fol- 
lowed by Co@585. The last one has a negative CP be- 
cause localized levels above the Fermi-energy scatter ma- 
jority spins, generating a current with excess of minority 
spin electrons. One should note that, despite the high 
magnetization of all Mn complexes, the absolute values 
of CP are not so remarkable in such systems. In partic- 
ular, the Mn@MV has a negligible CP, even though it 
has m = 3.0 fiB- For the noble metals, an unusual local 
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TABLE I: Binding energies (Eb) in eV, magnetization m in /is, and metal-carbon bound lengths (dmc) in A for all geometries 
considered in this work. (Imc are shown between brackets when there are different lengths. When convenient, we show either 
the position (see Fig. 1) or the local symmetry group. The binding energy is defined as Ef"'"'^'^''^'"'^ = E^^"^'^'^''^'"'^ - Ef^°^ - 

■^defect 
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magnetic moment generated by the polarization of "s" 
atomic orbitals leads to CP smaller than 10%. These 
results indicate that the existence of a local magnetic 
moment is a necessary but not sufficient condition to 
have a large CP. This can be clearly seen in the Fig. 
[3] (b) , where we show the CP as a function of the local 
magnetic moment, showing that there is no correlation 
between the local magnetic moment and the current po- 
larization. In Fig. [3] (c) we show the CP as a function 
of the magnetic moment induced only in the graphene 
carbon atoms. This magnetization occurs due to the hy- 
bridization between the graphene states and the adatom 
orbitals. However, this hybridization may occur (i) at 
deep levels, not affecting the CP, or (ii) near the Fermi 
energy, affecting the transmission probability near the 
transport window. Therefore we conclude that the CP 
behavior is not governed solely by the magnetic moment 
induced in graphene, but depends fundamentally on the 
energy window where the hybridization responsible for 
the appearance of such a magnetic moment occurs. 

In order to understand the ruling factors that govern 
the current polarization, we performed band structure 
calculations. From these results, we conclude that the 
degree of current spin polarization is governed by (i) 
the presence of localized levels, non degenerated in spin, 
nearby the Fermi level and (ii) by their hybridization 
with graphene levels. To illustrate these conclusions we 
present in detail three representative cases shown in Fig. 
m In panel (a) we plot the energy bands for the Au@MV 
complex. In this situation, there is a flat band around 
0.2 eV below the Fermi energy, with majority spin asso- 
ciated with Au "s" atomic orbitals, which is responsible 
for the local magnetic moment of 1.0 fiB- This band 
does not modify the shape of the Dirac cone because it 
has a weak coupling to the graphene tt bands. As a con- 
sequence, the transmission probability [T{E)] shown in 
Fig. |4] (b) is very similar to the pristine case around 
the Fermi energy. This example shows that a magnetic 
moment and localized states near the Fermi energy are 
a necessary, but not a sufficient, condition to generate 



spin-polarized currents in graphene, as discussed above. 
In fact, to generate current polarizations it is also neces- 
sary to have a high coupling between the localized levels 
and graphene it bands, modifying in this way the trans- 
mission probability for only one or both spin channels. 
As an example of modifications in only one spin channel, 
we show the energy bands for the Co adsorbed on pristine 
graphene in Fig. |4] (c) . It is notable that a Dirac cone is 
well defined only for the majority spin. For the minor- 
ity spin, localized levels hybridize with the Dirac cone 
states, and hence significantly modify the dispersion of 
these states. As a result, the transmittance [shown in 
Fig. 2] (d)] is very similar to the one in pristine graphene 
only for the majority spin, causing in this way a high cur- 
rent polarization. For the iVi@555777 complex there are 
modifications in the band structure for both spin chan- 
nels [see Fig. 0] (e)]. In this case, the local magnetic 
moment comes from localized states around the Fermi 
energy with a high hybridization with the Dirac cone. 
As a consequence, the transmittance for both spin chan- 
nels, shown in Fig. 111(f), are not only distinct from the 
pristine case, but also very different between them. In 
this system we still observe a flat band at 0.75 eV above 
the Fermi energy. However, as there is a weak coupling 
with the Dirac cone states, there is almost no change 
in T{E). Therefore, the examples shown in Fig. 2] elu- 
cidate the conditions to have spin-polarized currents in 
graphene generated by doping with magnetic impurities. 



C. Dependence on the doping concentration 

In most calculations we performed the distance be- 
tween the metallic adatom and its periodic image 
is 13.0 A. In order to understand the effect of the 
doping concentration we calculated the transmittance 
for Mn@MV also for supercehs of 17.3, 21.6, and 
26.0 A wide, as shown in Fig. [SJ The main charac- 
teristic of Mn@MV transmittance is a valley in the mi- 
nority spin transmittance around 0.5 eV above the Fermi 
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FIG. 4: (Color online) Band structure and spin-resolved 
transmittance for (a) and (b) Au'^MV , (c) and (d) Co over 
pristine graphene, (e) and (f) Ni@555777. The solid (dashed) 
lines correspond to the majority (minority) spin channels. For 
the sake of comparison, we also show the transmittance for 
pristine graphene. The smooth behavior of all curves indi- 
cates the convergence, regarding the fcx-points, to the true 
2D-nature of graphene. In the energy bands, we also indi- 
cate the position of the Dirac cone (with a vertical dashed 
line) , corresponding to the K symmetry point in a hexagonal 
supercell. 



energy, and another valley for the majority spin transmit- 
tance around 0.7 eV below the Fermi energy. With the 
widening of the supercell, the valleys became more local- 
ized, as a reflect of the decrease in the energy dispersion 
of the defect levels. However, the main results obtained 
with the smaller supercell arc unchanged. 
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FIG. 5: (Color online) Spin-resolved transmittance for 
MnQMV. The distance between Mn atoms are 13.0, 17.3, 
21.6 and 26.0 A in (a)-(d), respectively. 



potential the current polarization is negligible. However, 
with the inclusion of a negative gate the current has ex- 
cess of majority spin electrons, whereas with a positive 
gate there are excess of minority spin electrons. This 
behavior occurs because the gate potential shifts the lo- 
calized states that hybridize with the graphene tt bands, 
as shown in Figs. [nijb)-|ni^f), allowing a tuning of the val- 
leys in the transmittance (indicated by vertical lines) . In 
Fig. Elg), we also present CP as a function of Vgate- It is 
remarkable that CP can be changed from 50% to —18% 
with the application of a gate. Therefore, without any 
external magnetic perturbation, it is possible to control 
the degree of spin polarization in such systems. 
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D. Effect of a gate potential 

We also show that it is possible to control the current 
polarization via the effect of a gate potential. As an ex- 
ample, we show in Fig. |6]the transport properties for the 
MrMMV complex, as a function of Vgate- In panel (a) 
we present the spin-resolved currents. Without a gate 



FIG. 6: (Color online) (a) Spin-resolved current for the 
Mn@MV complex, were the circles (squares) represent the 
majority (minority) spin, (b)-(f) Transmission probability for 
Vgate = 3.0, 1.0, 0.0, —1.0 aud —3.0 V , respectively. Vertical 
lines indicate the position of the valleys in the transmittance. 
The Fermi energy is fixed by the infinite graphene leads, (g) 
CP as a function of Vgate- 
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IV. CONCLUSION 

In conclusion, via studies of adsorbed metal atoms in 
graphene with and without defects, we show that the 
presence of a local magnetic moment is a necessary but 
not sufficient condition to create a current polarization. 
Essentially, a fundamental requirement to have a non- 
zero current polarization is the presence of spin-split lo- 
calized levels near the Fermi energy that strongly hy- 
bridize with the graphene Dirac-cone states. We also 



show that with an external gate potential it is possible 
to control the degree of current polarization without any 
magnetic perturbation. The knowledge of such proper- 
ties for isolated scatterers is a fundamental first step to 
design spintronic devices, where other important effects 
such as multiple scattering and disorder may even en- 
hance the polarization signal^"—. 
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